Activated escape of periodically modulated systems 
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The rate of noise- induced escape from a metastable state of a periodically modulated overdamped 
system is found for an arbitrary modulation amplitude A. The instantaneous escape rate displays 
peaks that vary with the modulation from Gaussian to strongly asymmetric. The prefactor v in the 
period-averaged escape rate depends on A nonmonotonically. Near the bifurcation amplitude A c it 
scales as v oc (A c — Ay. We identify three scaling regimes, with £ = 1/4, — 1, and 1/2. 
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Thermally activated escape from a metastable state is 
often investigated in systems driven by time-dependent 
fields. Recent examples are activated transitions in mod- 
ulated nanomagnets Q, 0, and Josephson junctions 
0,0,0- Modulation changes the activation barrier. This 
enables efficient control of the escape rate and accurate 
measurement of the system parameters jjj |g| . 

Most frequently used types of modulation are slow 
ramping of a control parameter, when the system re- 
mains quasistationary, and periodic modulation. In the 
latter case the system is away from thermal equilibrium, 
which complicates the theoretical formulation of the es- 
cape problem . 

In the present paper we extend to periodically mod- 
ulated systems the analysis of the escape rate done by 
Kramers for systems in thermal equilibrium [10j. Our 
approach gives the full time-dependent escape rate W(t) 
as well as the period-averaged rate W = ^exp(— R/D), 
where R is the activation energy of escape and D is the 
noise intensity, D = ksT for thermal noise. 

For comparatively small modulation amplitude A es- 
cape of an overdamped Brownian particle was studied in 
Ref.lHl The range of intermediate A and intermediate 
modulation frequencies ujf was analyzed in Refs. Il2l IT^L 
Here we find W(t) for an arbitrary A and an arbitrary 
interrelation between up and the relaxation time of the 
system t T . We show that the prefactor v depends on A 
strongly and nonmonotonically. It displays scaling be- 
havior near the bifurcational modulation amplitude A c 
for which the metastable state disappears. 

In the spirit of Kramers' approach, we relate the in- 
stantaneous escape rate W(t) to the current well behind 
the boundary (?&(<) of the basin of attraction to the ini- 
tially occupied metastable state (q is the system coor- 
dinate). This is the current usually studied in experi- 
ments. Because of the oscillations of qb(t), it has a differ- 
ent functional form from the current at qb(t) calculated 
in Refs. EUll 

We find W(t) by matching the probability distribu- 
tion p(q,t) near qb{t) and inside the basin of attraction. 
This can be done in a broad parameter range without a 
complete calculation of p(q,t) near cfa(t), using singular 
features of the dynamics of large fluctuations. 

For a periodically modulated overdamped Brownian 
particle, the distribution p(q, t) is given by the Fokker- 



Planck equation (FPE) 



d t p=-d q [K(q,t)p\+Dd 2 qP . 



(1) 



Here, K(q,t) is the periodic force driving the particle, 
K(q, t) = K{q, t + tf) = —d q U{q 1 t), where tf — 2-k/ujf 
is the modulation period and U (q, t) is the metastable 
potential. The equation of motion of the particle in the 
absence of noise is q = K(q,t). The metastable state 
q a (t), from the vicinity of which the system escapes due 
to noise, and the basin boundary qb(t) are the stable and 
unstable periodic solutions of this equation, respectively. 

We will assume that the noise intensity D is small. 
Then in a broad time range t r <C t <C X/W the distribu- 
tion p(q, t) is nearly periodic in the basin of attraction 
to q a (t). The current away from this basin, and thus the 
escape rate W{t), are also periodic. 

The distribution p is maximal at q a (t) and falls off 
exponentially away from it. In the presence of periodic 
driving it acquires singular features as D — > UM . some 
of which have counterparts in wave fields |15| . with D 
playing the role of the wavelength. The singularities ac- 
cumulate near (ft,(t). In order to find W(t) one has to 
understand how they are smeared by diffusion. 

In the absence of noise the motion of the system close 
to the periodic states qi(t) (i = a,b) is described by the 
equation q = K with K linearized in q — qi(t). The 
evolution of q(t) — qi(t) is given by the factors 



Ki(t, t') = exp 



drpiir) 



(i = a,b), 



(2) 



where Pi(t) = pi(t + tf) = [d q K(q,t)] qi ny Over the 
period tf the distance q(t) — qi{t) decreases (for i = a) 
or increases (for i = b) by the Floquet multiplier Mi = 
Ki(t + Tp,t) = exp(p,iTF), where pt is the period-average 
value of Pi(t), with p a < 0, pb > 0. 

For weak noise the expansions of K can be used to 
find p(q,t) near q a ,b(t). Near the metastable state q a , 
the distribution is Gaussian [l6| . p(q,t) oc exp{ — [q — 
q a (t)] 2 /2Da^(t)}. The reduced time-periodic variance is 
given by the equation 

a l 2 (t) = 2|Mr 2 -ip 1 f * dt lK - 2 {t + t u t) (3) 
Jo 

with i = a (in the absence of modulation a 2 — l/\p a \). 
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The general form of the periodic distribution near the 
unstable state %{t) (the boundary-layer distribution) can 
be found from Eq. using the Laplace transform, simi- 
lar to the weak-driving limit With K linear in q — q b , 
the equation for the Laplace transform of p{q, t) is of the 
first order, giving 

/>OC 

p(q,t) = / dpe^' D p{p,t), Q = q-q b (t), 
Jo 

p(p,t) = ED- 1 ' 2 exp{- [a(4>) + p 2 o 2 (t)/2] /D) .(4) 

In Eq. £ is a constant, s(</>) is an arbitrary zero-mean 
periodic function, s(</> + 2tt) = s(<f>), and <f> = 4>{p, t), 

0(p,t) - n F ln\pK b (t,t')/fi b l D ]. (5) 

Here, Sip = ujp/p b = 27r/lnMb is the reduced field fre- 
quency, Id = (2D/p, b ) 1 ^ 2 is the typical diffusion length, 
and t' determines the initial value of (/>; from Eq. JSJ), 
<j>(p,t + Tp) — <fi(p,t) + 2tt. In Eq. Q we assumed 
that the basin of attraction to q a lies for q < q b (t), and 
\Q\ < Aq = mint [q b (t )-?«,(*)]• 

The form |0J is advantageous as it immediately gives 
the current j(q, t) from the occupied region (— oc, q\. Well 
behind the basin boundary, where Q = q — q b (t) ^> Id, 
diffusion can be disregarded, the current becomes con- 
vective and gives the instantaneous escape rate. With 
q « q b + p b Q, we have j(q, t) w p b (t)p(q, t)Q, at a given 
Q. Disregarding the term oc p 2 /D in p for Q 3> Id, we 
obtain from Eq. Q 

/>oo 

j(q,t) = fib^ED 1 / 2 dxe- x exp[-s{<f> d )/D]. (6) 
Jo 

Here, 0d = Qf^[x K b (t d ,t')], and t d = t d (Q,t) is given 
by the equation n b (td,t) = Id/2Q. In the whole har- 
monic range j depends on the observation point Q only 
in terms of the delay time td, which shows how long it 
took the system to roll down to the point Q, dtd/dQ = 
— l/p b (td)Q- We note that p b (t) can be negative for a 
part of the period, leading to reversals of the instanta- 
neous current. 

The escape rate W is given by the period-averaged 
j(q,t) and is independent of q. From Eq. 

W=^SD 1 / 2 d0exp[-s(0)/£>]. (7) 
2n Jo 

Eqs. © and (JZJ) provide a complete solution of the 
Kramers problem of escape of a modulated system and 
reduce it to finding the function s. They are similar in 
form to the expressions for the instantaneous and average 
escape rates for comparatively weak modulation, \s\ <C 
R, where s was obtained explicitly 

Unless the modulation is very weak or has a high fre- 
quency, for small noise intensity max s ~ | min s\ 3> D. 
In this case the major contribution to the integrals in 
Eqs. ©, J7J comes from the range where s is close to its 
minimum s m reached for some (p — 0m- Then the escape 



rate j(q,t) sharply peaks as function of time once per 
period when 4>d(t) — 4> m - This means that escape events 
are strongly synchronized. As we show, both j(q,t) and 
W are determined not by the global shape of s(<fi), but 
only by the curvature of s(<p) near </> m . 

To find j(q,t) we match Eq. (0J to the distribution 
p{q,t) close to q b (t) but well inside the attraction basin, 
—Q 3> Id- For small D this distribution can be found, 
for example, by solving the FPE Q in the eikonal ap- 
proximation, p{q,t) — exp[— S{q, t)/D\. To zeroth or- 
der in D, the equation for S — So has the form of the 
Hamilton- Jacobi equation dtSo = —H for an auxiliary 
nondissipative system with Hamiltonian |17| 

H(q,p;t) = p 2 + P K{q,t), p = d q S . (8) 

The Hamiltonian trajectories of interest for escape 
q(t),p(t) start in the vicinity of the metastable state. The 
initial conditions follow from the Gaussian form of p(q, t) 
near q a (t), with S = [q - q a (t)} 2 /2a 2 a (t). 

To logarithmic accuracy, the escape rate is determined 
by the probability to reach the basin boundary q b {t), i.e., 
by the action So(q b (t), @- The Hamiltonian trajectory 
Qopt(t) , Popt(t) , which minimizes So(q b (t), i), approaches 
q b (t) asymptotically as t — > oo. This is a heteroclinic tra- 
jectory of the auxiliary system, it is periodically repeated 
in time with period tf', Q op t(t) is the most probable es- 
cape path (MPEP) of the original system. 

Close to q b (t), the Hamiltonian equations for q(t),p(t) 
can be linearized and solved. On the MPEP 

Popt(t) = -Qo P t{t)/al(t) = K^(t,t')p opt (t'), (9) 
S (q op t(t),t) = R - Q 2 opt (t)/2<j 2 (t), 

where Q op t(t) = q op t(t) - q b (t). The quantity R = 
•So (<fopt (i) > t) t ->oo * s tne act i vat i° n energy of escape. 

The surface So (q, t) is flat for small Q — Q op t due to 
nonintegrability of the dynamics with Hamiltonian © 
[Til . It touches the surface S b {q,t) = R - Q 2 /2a 2 {t) 
on the MPEP, Q = Q opt (t). Away from the MPEP 
So(q,t) > S b (q,t), and therefore the function p b (q,t) — 
p(q,t) exp[5f,(q, t)/D] is maximal on the MPEP. 

We match on the MPEP p b found in the eikonal ap- 
proximation to the maximum of p b found from Eq. Q 
near the basin boundary. For |s m | D and — Q ^ Id, 
the integral over p in Eq. Q can be evaluated by the 
steepest descent method. The integrand is maximal if 
p = —Q/o- 2 (t) and s is minimal for this p, i.e., <p(p,t) — 
(j) m and s = s m . These conditions can be met on the 
whole MPEP at once, because 4>(p op t(t), t) =const. Then 
from Eq. (gj 

p(q,t) =£ b {t)exp[-S b (q,t)/D], 

£ b (t) = SD- 1 ' 2 [a 2 (t) + Q 2 F OoptW] ~ V2 . (10) 

where £ = £(2ttD)^ 2 exp[(i? - s m )/D], and s m = 
[d 2 s / d(j) 2 ]^ m - From Eqs. 10, ijTU|) . not only the expo- 
nents, but also their slopes coincide along the MPEP for 
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the boundary-layer and eikonal-approximation distribu- 
tions. 

The function £b{t) should match on the MPEP the 
prefactor of the eikonal-approximation distribution p = 
exp(— S/D), which is given by the term S± oc D in S. On 
the MPEP, z = exp(25i/D) obeys the equation [Tllisl l 



d 2 z/dt 2 - 2d{zd q K)/dt + 2zpd 2 q K = 0, 



(11) 



where q = q op t{t), p = p op t(t). The initial condition to 
this equation follows from p(q,t) — z^ 1 / 2 exp(— Sq/D) 
being Gaussian near q a (t), which gives z(t) — > 2irDa 2 {t) 
for t — > — oo. Close to qb(t), from Eq. z(t) = 

D[Zia 2 (t) + Z 2 Po p 2 t (£)], where Zi :2 are constants [l2j : 
the term cx Z\ was disregarded in the analysis 0. Re- 
markably, z~ 1 / 2 (f) is of the same functional form near 
qb{t) as £b{t). Thus the prefactors in p(q,t) as given by 
the eikonal and the boundary-layer approximations also 
match each other. 

Explicit expressions for the escape rate in the regime of 
strong synchronization can be obtained for comparatively 



weak or slow modulation, where s ' 



Q, 2 F C < R. 



s m > D but 



(12) 



The results for D -C \s m \ <C R should coincide with the 
results of Ref. ITTL which were obtained in a different way. 
We have verified this by finding s'^ from Eq. 1(11(1 by 
perturbation theory in the modulation amplitude A. 

Condition 1(12(1 can be met for large A, where ~ R, 
provided the modulation frequency is small, uj F t r ~ 
tip <C 1 (adiabatic modulation). Here, the MPEP is 
given by the equation q opt — —K{q opt , t m j , with t m found 
from the condition of the minimum of the adiabatic bar- 
rier height AU{t) = U{q b {t),t) - U{q a {t),t). The acti- 
vation energy R — AU m = AU(t m ). 

The value of can be obtained from z(t) or 
by matching the adiabatic intrawell distribution oc 
exp[— U(q,t)/D] and the boundary layer distribution Q 
in the region \Q\ » lp> and tips'^ <C Pb{t m )Q 2 for 
I* — tm\ -C rjr. Both approaches give Clpfj^s'^ — AU m , 
where pb and AC/ m = dfAU are calculated for < = t m . 

The form oij{q, t) depends on the parameter tlps'^/D. 
When it is small, the term oc p~ pt in £b(t) [Eq. 1(10(1 ] and 
z(t) is also small away from the diffusion region around 
qb. Then z = 2irDa 2 {t m ). The pulses of j(q,t) are Gaus- 
sian, 



I Ma Mb 



1/2 



2?r 



e -fl./Z5V^ e -(t~t k ) 2 AU m /2D 



(13) 



[Ma,6 = Ma,h(*m)]- They are centered at i fe = £ m + 
fcrF, with fc = 0, ±1,... [we disregard the delay ~ 
p^ 1 h\{Q /Ip,) in tk}- Eq. 113|) corresponds to the fully 
adiabatic picture, where the escape rate is given by the 
instantaneous barrier height AU(t). 

The current has a different form for tips'^/D ^> 1. 
Because p~ 2 t (t) oc Kj(t, t') exponentially increases in time 
near qb, the term oc p~ 2 t in £b and z becomes dominating 



before the MPEP reaches the diffusion region \Q\ ~ 
Then Eqs. (fTU|) give 



£fc = exp(27rfc/r2F), 



ar fc e- a *, (14) 



k— — oo 

so = p op t(t)Q/D 



Note that here p pt (i) can be smaller than lp> / a 2 (t) . 

Eq. 1(14(1 describes the escape rate in the whole region 
^F s m ^ -^i it does not require the adiabatic approx- 
imation. Its form is totally different from that of the 
diffusion current —Ddqp on the basin boundary Q — 

as given by Eqs. @, (TO). The ratio = £l F Z~ 1/2 

can be obtained by solving Eq. (|lip. 

For ili? <C 1 the current (|14ll is a series of dis- 
tinct strongly asymmetric peaks, with Xk ~ exp[— (t — 
krp — t m )fib(tm)] near the maximum. The transi- 
tion between the pulse shapes 1(13(1 and 1)14(1 occurs for 
tips'^/D ~ 1. It is described by Eq. with E = 
(2TT)- 1 D- 1 / 2 \p a /p b \ 1/2 exp[-(R-s m )/D}. We note that, 
for flp <C 1, the shape of current pulses in the whole 
range (O is the same as for weak modulation |l9j . but 
the parameters depend on A,lof differently. 

With increasing Hp the peaks of j (|14l) are smeared out 
and the escape synchronization is weakened. For flp 3> 1 
it disappears (s^ rapidly decreases with u)p for large £l F ). 

The escape current ((1411 is completely different from 
the current on the basin boundary 12]. The regime 
rt F s'^/D <C 1, where the current has the form l|13|l . can- 
not be studied in the approximation [l2T | at all. 

In the range ~ |s m | ^> D, the period-averaged es- 
cape rate (0) is 

W = vexp(-R/D), v = p b £D 1/2 /2TTy^. (15) 



The prefactor v can be expressed in terms of Zo , formally 
giving the result 01 even where the theory ^3] does not 
apply. 

The asymptotic technique developed in this paper al- 
lows obtaining the prefactor v in several limiting cases. 
For comparatively weak modulation, D <C | s m j <C R, 
Eqs. (fTTJl . itTB^I give the same result as in Ref. Since 
the theory [ll| covers the whole range \s m \ <C R, a transi- 
tion from the Kramers limit of no modulation to the case 
of arbitrarily strong modulation is now fully described. 

In the whole range where the adiabatic approximation 
applies, fl F <C 1, we obtain 



v = (2n)-^ 2 \p a p b \ 1 / 2 D 1 / 2 LUp(AU nl )- 1/2 



(16) 



where p a .b are calculated for t — t m . Interestingly, v i|ltj|) 
is independent of the modulation frequency. 

Eq. ((16(1 is simplified for the modulation amplitude 
close to the bifurcational value A c where the metastablc 
and unstable states merge. For small SA = A c — A 
and wp\t — t m \ the adiabatic barrier is AU(t) oc [SA + 
a c u 2 F (t - ■U i ff /2 (hereo, ~ A c ), and |yu aj6 | oc {SA) 1 / 2 , 
cf. Refs. Hfl, HH H2, El, HI, HE Then, from Eq. TO, the 
prefactor in the adiabatic limit scales as v oc {SA) 1 / 4 . 



4 




0.15 A 0.3 



FIG. 1: The prefactor v in the average escape rate W CHJ- 
The results refer to a Brownian particle with K(q, i) — 
q 2 — 1/4 + Acos(u>Ft), lof = 0.1, and describe escape in the 
regime of strong synchronization, where v oc D 1/2 . The solid 
line for small A show the scaling v oc A~ 1//2 11]. The solid 
line for small A c — A (A c w 0.29) shows the scaling 1171 . 
The dashed line shows the result of the numerical solution of 
Eq. 1111 . The squares and crosses show the results of Monte 
Carlo simulations for R/D — 5 and R/D — 6, respectively. 

The slowing down of the system motion makes the adi- 
abatic approximation inapplicable in the region SA/A C < 
ftp. In contrast to the adiabatic scaling R qc_(cL4) 3 / 2 , the 
activation energy scales here as R oc {SA) 2 24] . Using the 
results [24|, we obtain from Eq. (|llfl 

v= V0 D 1 '\SA)- 1 u*J*, (17) 

where v = (G^Wp)- 1 /^^ d^K] 1 /* /\d A K\. Here all 
derivatives are evaluated for q, t, and the amplitude A = 
A£ d A c where the minimum and maximum over q of 



the potential U(q,t) merge (once per period). 

From Eq. I|17|) . the prefactor v oc (SA)^ 1 sharply in- 
creases as the modulation amplitude approaches A c . This 
is qualitatively different from the decrease of v in the 
adiabatic approximation. The scaling v oc (SA) -1 agrees 
with the numerical solution of Eqs. I|llfl . 115fl for a model 
system shown in Fig.^ The calculations in a broad range 
of A are also confirmed by Monte Carlo simulations. 

For high frequencies, f2f 3> 1, modulation does not 
lead to exponentially strong escape synchronization. The 
prefactor in the escape rate is v — |/Z„/i&| 1 ' 2 /27r, it is in- 
dependent of the noise intensity D. Near the bifurcation 
point it scales as in stationary systems, where v oc (SA) 1 / 2 
and R oc (SA) 3 / 2 [U EH Very close to the bifurcation 
point modulation is necessarily fast, because \p, a .b\ —> 
for A — > A c . Therefore the prefactor always goes to zero 
for A — > A c . However, for small ujp the corresponding 
region of SA is exponentially narrow p4| . 

In conclusion, we have obtained a general solution of 
the problem of noise-induced escape in periodically mod- 
ulated overdamped systems. With increasing modulation 
frequency, the pulses of escape current change from Gaus- 
sian to strongly asymmetric; for large up current mod- 
ulation is smeared out. The prefactor v in the period- 
averaged escape rate is a strongly nonmonotonic func- 
tion of the modulation amplitude A for low frequencies. 
It first drops with increasing A to v oc (D/A) 1 / 2 [Tl| . 
then varies with A smoothly [uH H^ ] , and then sharply 
increases, v oc D 1 / 2 / (A c — A) near the bifurcation ampli- 
tude A c . We found three scaling regimes near A c , where 
v oc (A c - A)( with C = 1/4, -1, or 1/2. The widths of 
the corresponding scaling ranges strongly depend on the 
modulation frequency. 

This research was supported in part by the NSF DMR- 
0305746. 
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